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ABSTRACT 

In a previous paper we described a method of estimating the single-measurement bias to be 
expected in astrometric observations of targets in crowded fields with the future Space Interfer- 
ometry Mission (SIM). That study was based on a simplified model of the instrument and the 
measurement process involving a single-pixel focal plane detector, an idealized spectrometer, and 
continuous sampling of the fringes during the delay scanning. In this paper we elaborate on this 
"instrument model" to include the following additional complications: spectral dispersion of the 
light with a thin prism, which turns the instrument camera into an objective prism spectrograph; 
a multiple-pixel detector in the camera focal plane; and, binning of the fringe signal during scan- 
ning of the delay. The results obtained with this improved model differ in small but systematic 
ways from those obtained with the earlier simplified model. We conclude that it is the pixellation 
of the dispersed fringes on the focal plane detector which is responsible for the differences. The 
improved instrument model described here suggests additional ways of reducing certain kinds of 
confusion, and provides a better basis for the evaluation of instrumental effects in the future. 

Subject headings: Astronomical Instrumentation 



1. Introduction 

The Space Interferometry Mission PlanetQuest 
(hereafter SIM) is being designed by NASA/JPL 
to carry out a program of extremely precise as- 
trometry on stars, thereby contributing to a wide 
variety of research topics in astronomy includ- 
ing the study of the mass distribution in the 
Galaxy and the search for earth-like extra-solar 
planets. A few of these topics involve making 
precise position measurements of target stars in 
areas of the sky containing more than w 0.4 
detectable stars per square arcsecond. In these 
"crowded-field" cases we have previously shown 
( Sridharan fc AllenI 2007 . hereafter Paper I) that 
the individual astrometric measurements made on 
such targets could be biased by the light from ex- 
traneous field stars. In Paper I of this series, we 
introduced a simplified instrument model for SIM, 
and used it to evaluate the likelihood of confusion 
bias on a number of Key Projects to be carried 



out as part of the SIM science program. In that 
paper we concluded that, in the small number of 
cases where confusion bias may be problematic, 
the likelihood of bias can be reduced by judicious 
design of the observing program. 

The measurement model used in Paper I is 
based on a phasor description for estimating the 
contributions of all stars in the field of view 
(FOV). This part of our model is not in question 
here. However, there were three subsequent sim- 
plifications to the instrument model which were 
suspect: First, the detector in the focal plane was 
assumed to consist of one large pixel per chan- 
nel which collected all the light diffracted through 
the field-stop; second, the details of just how the 
light is dispersed in wavelength were ignored; and 
third, the total fringe signal was assumed to be 
continuously sampled during the scanning of the 
delay. In this paper we elaborate further on the 
instrument model to include the following com- 
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plications: A multiple-pixel detector in the focal 
plane; spectral dispersion of the light with a thin 
prism, which effectively turns the instrument into 
an objective prism spectrograph; and, integration 
and discrete sampling of the fringe signal during 
scanning of the delay. We then repeat the com- 
putations of confusion bias for a few specific cases 
already evaluated with the simplified approach 
of Paper I. In general, the improvements to the 
instrument model introduced in this paper have 
only a minor impact on the biases, and come at 
the cost of significant complication. However, we 
have not yet implemented a number of instrumen- 
tal effects which may be present in real SIM data, 
including pointing errors, delay line jitter, channel 
band-shape changes, and so on. Besides its ped- 
agogical value, the more complete description of 
how SIM works which we give here will be useful 
for such studies in the future. 

We begin with stepwise description of the oper- 
ation of SIM, progressively adding complexity at 
each step. This leads us to the concept of the fan 
diagram, which was introduced in the context of 
SIM in Paper I as an aid to understanding how 
SIM works, and which will gain further in impor- 
tance here. We then use this improved instrument 
model to re-evaluate the confusion bias for several 
cases presented in Paper I. Finally, we point out 
a new possibility related to the pixellation of the 
focal plane which may be used to further reduce 
the likelihood of confusion bias from field stars lo- 
cated within the SIM FOV but separated by > 1" 
from the target star. 

2. SIM Beam combination, dispersion, 
and focal plane imaging 

We start by considering how the light is pro- 
cessed after the beams from each arm of the in- 
terferometer are combined. As described in Ap- 
pendix A of Paper I (see especially Figure 19), the 
collimated beams from each collector are passed 
through the internal delay lines, and the two 
beams are superposed in the beam combiner. The 
precise operation of these optical assemblies is not 
directly relevant to our present goal, so we will 
not elaborate further on them here. The combined 
beams then pass through the prism disperser, and 
are imaged on a fast-readout charge- coupled de- 
vice (CCD) in the focal plane of the camera. Fig- 



ure [T] is a sketch of this optical system. 

In order to understand the structure of the im- 
ages which appear in SIM's focal plane, it is use- 
ful to proceed in a series of steps of increasing 
complexity. For the moment, we will assume that 
SIM's primary apertures are very large, so that 
star images appear as unresolved dots (this will 
be rectified later). Figure [2] shows a sketch of the 
2-D images which appear in such an "ideal" focal 
plane (at the location of the detector at the right 
side of Figure [1]). 

• If one of the siderostats is covered and the 
dispersing prism removed from the light 
beam, the image appearing at the focal plane 
will show the stars in the FOV at their re- 
spective positions. In the top panel of Figure 
[2]we have sketched a field with 3 stars of dif- 
ferent spectral types in the FOV: # 1 is a G 
star at the center; # 2 an O star to the upper 
left; and, # 3 an M star to the lower right. 

• Next, if we introduce the prism into the light 
path, the images of the stars will be "spread 
out" in the direction of the dispersion (as- 
sumed to be the +X direction in Figure[T]) in 
the same fashion as with an objective prism 
spectrograph; the focal plane will consist of 
"patches" of light, each varying according to 
the SED of its star, as shown schematically 
in the middle panel of Figure 

• If we now uncover the second siderostat, 
each star's spectrum will become an inter- 
ferogram consisting of its SED modulated 
by a complex pattern of fringes depending 
on the setting of the internal delay and the 
location of the star on the sky. This modu- 
lation pattern can be read off along vertical 
lines in the fan diagram of Figure [3] (to be 
described in the next section). 

• Finally, if we scan the internal delay (5, each 
star's interferogram will change in a different 
way, according to location of the star on the 
X-axis of the fan diagram. If the white fringe 
happens to be positioned exactly at the lo- 
cation of a star, the interferogram of that 



^The SEDs are shown here as graphs for simphcity; in fact 
they will appear as elongated patches of light with intensity 
distributions in +X as shown in the graphs. 
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Fig. 1. — Schematic of the optics that follow the beam combiner and which are responsible for dispersing and 
imaging the combined light from the two collectors. The combined beams from the two siderostat collectors 
enter at the left and proceed through the compressor and field stop to the prism disperser, and are finally 
imaged on the detector at right. 



star will lose its fringes and show simply the 
stellar SED. 

As an aid to understanding the final point 
spread function for SIM which results from the 
steps described above, we recall the concept of the 
fan diagram which was introduced in Paper I. 

3. The Fan diagram 

The response pattern P{S,d,k) of SIM's as- 
trometric (Michelson) interferometer in a quasi- 
monochromatic channel can be written as (cf. Pa- 
per I): 

P{S,e,k) = Pa{l + Asin{27rk[6 ^ B9])}, (1) 

where Pq is the total power, A is the resultant 
fringe amplitude, k is mean wavenumber of the 
channel, S is the internal path delay, and B9 is the 
external path delay where B is the projected base- 
line length and 6 is an angle on the sky measured 
from the direction perpendicular to the projection 
of the interferometer baseline. A contour plot of 
this response pattern as a function of the optical 
path difference (OPD) A = 5 - 56* and A = 1/fc 
is shown in Figure [3l It has the appearance of a 
hand-held collapsable fan, and hence the name fan 
diagran^ The fan diagram is a convenient and 

^In their pape r on observing binary stars with SIM, 
iDalal &: Griestj ll200ll ) introduce the "channeled spectrum" , 
which is related to the fan diagram defined here. 



compact way of displaying the response of SIM's 
interferometer in terms of delay and wavelength. 
Although it is not itself a "picture" of SIM's focal 
plane, as we shall see it is an aid in constructing 
such a picture for an arbitrary ensemble of target 
and field stars. A few salient features of the fan 
diagram are: 

1. It is an intrinsic property of the instrument; 
it describes the theoretical point spread func- 
tion of SIM, i.e. the idealized response to 
a point source at position 9 at mean wave- 
number k and setting 6 of the internal delay. 

2. A horizontal profile through the diagram 
shows a single fringe with a period which 
increases gradually from short to long wave- 
lengths. 

3. The diagram shows an (anti-)symmetry 
about the zero point on the x axis. At this 
OPD, the fringe position is the same at all 
wavelengths, hence this position is called the 
white fringe. At larger x offsets, the fringes 
are increasingly tilted, and the "blades of 
the fan" become more and more horizontal. 

4. A star in the FOV is represented by a verti- 
cal line in the fan diagram. We may think of 
this vertical line as shifting along the x axis 
as the internal delay S is changed. The ba- 
sic astrometric measurement procedure on a 
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Fig. 2. — Schematic of the "ideal" focal plane. Top: When the light from one of the siderostats of SIM is 
blocked and the prism removed from the light beam, the stars will appear as spots of light in this idealization. 
Stars 2 and 3 are offset from star 1 in this sketch by w 1" in each dimension. Middle: The star images 
will spread into spectra when the prism is introduced into the path of the (single) beam. The three stars 
1, 2, and 3 shown for illustration in the field arc of spectral type G, O, and M, respectively; note how their 
SEDs differ. We have adopted the convention that the short- wavelength end of the spectrum appears to the 
left (small X), and the wavelength increases with position X to the right; this is otherwise arbitrary, and 
depends on the design of the instrument. For this sketch we have further assumed that the relation between 
X and A is simply linear. Bottom: When the beam from the second siderostat is added, the spectra will be 
modulated with a complicated fringe pattern, as described in the text. If the white fringe of the instrument 
happens to coincide with one of the stars (assumed to be star 1 here), the interferogram of that star will 
show no fringes, but will reproduce the spectral energy distribution of the star. Stars 3 and 2 are at delay 
offsets of ±43 /um from star 1, respectively. 



single target star involves adjusting the in- 
ternal delay 6 until the white fringe coincides 
with the target delay position, i.e. until the 
vertical line is shifted to A = 0. 

5. The presence of multiple stars in the FOV 
will lead to a complicated response as stars 
at different values of each contribute dif- 
ferent fringe patterns. The total response is 
the sum of all these patterns. 

4. Modeling the real focal plane 

The real focal plane of SIM differs from the ide- 
alized model we have just presented in a number 
of important ways: First, owing to diffraction by 



the finite size of the collector, the images of stars 
will be large and will grow (in size) with increasing 
wavelength. Second, the dispersion by the prism 
may result in overlap of the diffracted images of 
different stars in a crowded field, and these images 
will not overlap at the same wavelengths. Third, 
the pixels of the CCD detector are much larger 
than those of our ideal model. Finally, the outputs 
in each CCD pixel are binned as the internal delay 
is scanned in order to measure the fringe param- 
eters (total power, amplitude and phase) in each 
wavelength channel. In this section we consider 
each of these additional complications in turn. 

To begin, we shall need a numerical model of 
the ideal focal plane which will allow us to ade- 
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Fig. 3. — The "fan diagram" , a contour plot of SIM's astrometric interferometer response (the "point spread 
function" , or PSF) as a function of two variables: The optical path difference (OPD) A = (5 — B6 between 
the internal and the external path delays, and; the wavelength A. The fringes appear to 'fan out' at longer 
wavelengths, hence the name. The point of (anti-)symmetry on the x-axis is called the white fringe. The 
vertical line at A = 5 — B6 = +1 /im traces the interferometer's spectral response to a star located in the 
FOV at an angular offset of —0.023" from the white fringe. Scanning the internal delay 6 will result in 
moving the location of this vertical line along the x axis. An "astrometric delay" measurement with SIM 
begins by adjusting S until the star is situated at the white fringe. The dotted, solid and dashed contours 
represent contour levels of 0, 1 and 2 respectively. 



quately represent the various smoothing and over- 
laps described above. This model involves adopt- 
ing a 2-D grid of data points with adequate res- 
olution in each direction. The various consider- 
ations which go into this choice are discussed in 
Appendix |Al we have settled on an array of size 
2112 X 192, and taken each element to represent 
the brightness in an area of 1 x 1 /im in the focal 
plane, corresponding to a square of side 1/12 arc- 
sec on the sky for scale expected to be adopted in 
SIM's camera. 

4.1. Diffraction 

SIM's siderostat collectors will have effective 
outer and inner diameters of 304.5 mm and 
178 mm, respectively. The point spread function 
will therefore be approximately an Airy function 
with a full width at half maximum which increases 
from w 0.3" at 400 nm to w 0.7" at 1000 nm. The 



field of view (FOV) of the sky which appears on 
the detector is limited by the field stop (Figure 
[1]); the stop size is expected to be equivalent to 
a circle of diameter 3" in the camera focal plane. 
Our model includes light diffracted through this 
stop from all stars lying within a circle of diameter 
6". 

4.2. Dispersion 

Next we must specify the prism dispersion, i.e. 
the relation between A and X. For a star at the 
center of the FOV, the prism used in SIM has a 
measured dispersion as shown in the left panel of 
Figure ID This relation can be closely approxi- 
mated with a polynomial; the coefficients are given 
in Appendix B of Paper I. It is likely that these 
coefficients will be somewhat dependent on the lo- 
cation of the target star in the field of view; how- 
ever, these details are unknown at present, so this 
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Fig. 4. — Left: Experimentally-measured dispersion data for the current design of SIM's prism disperser. 
The dotted line is for wavelength (right axis) and the solid line for wavenumber (left axis). Note that the 
dispersion is a non-linear function of wavelength A, but close to a linear function of wavenumber k. Right: 
Estimated throughput of the entire instrument including reflectivity of the optics and the detector response. 
The dotted curve is a convenient polynomial fit to the measured values (solid curve). 
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Fig. 5. — Model of the "ideal" image in the focal plane of SIM for a single target star, showing the combined 
effects of diffraction and dispersion. Note that the size of the PSF increases as the wavelength increases 
in the +X direction (horizontally to the right). The intensity has been displayed on a logarithmic scale in 
order to show the "Airy rings" , which are no longer rings at all, but long horizontal features which spread 
out with increasing wavelength. For simplicity, a flat SED has been used in this illustration and hence the 
intensity variation in the +X direction mimics the throughput variation. 



effect has not been included in our models. 

The combined effects of diffraction and disper- 
sion in SIM's focal plane are shown in Figure \5\ 
Note that one result of this method of dispersing 
the fringes is that each point in the focal plane 
receives photons from a range of wavelengths of 
the target. The intensity at a given position in 
the focal plane will change, as will the nominal 
wavelength associated with that position, and the 
extent of this change will depend on the SED of 
the target. This will also lead to small changes 



in the estimated fringe phase. Note further that 
this effect is always present, even if there are no 
other stars within the FOV of SIM; it is a "self- 
confusion" effect. We will quantify this effect for 
SIM later in this paper. 

4.3. Throughput and fringe modulation 

There are two remaining modifications to the 
distribution of intensity in the focal plane. The 
first is caused by the variation in the overall 
throughput of the system as a function of wave- 
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length, including the reflectivity of the optics and 
the sensitivity of the CCD detector in the focal 
plane. The current best estimate of this through- 
put is shown in the right panel of Figure S) A 
polynomial fit to this function is given in Appendix 
B of Paper I. As with the wavelength dispersion, 
the throughput is also expected to be weakly field- 
dependent, but as yet we do not have sufficient 
knowledge to model this effect. 

The final (and in many ways most important) 
effect to include is the modulation of the spectra 
by the fringe modulation function. This is strongly 
field dependent. This function will differ for each 
star in the FOV, and is given by the fan diagram 
of Figure [31 which depends on the baseline length 
B, the location 9 of the star in the FOV, and the 
position of the internal delay d. 

A few clarifying remarks on the geometry are in 
order here. The relative angular location of each 
star in the FOV is measured as projected onto the 
interferometer baseline; see Figure 3 in Paper I 
for a sketch. The prism disperser in SIM's optical 
train is oriented such that the direction of disper- 
sion is parallel to the projection of the baseline on 
the sk30, and the camera CCD is oriented such 
that the +X pixel grid direction is also aligned 
with the baseline projection. Both of these lat- 
ter choices are a matter of convenience and not of 
necessity. We have used them here as well, since 
they simplify the model calculations. 

Figure [6] shows the ideal focal plane for two 
stars separated by « 1.5" and includes all effects 
discussed so far; diffraction, dispersion, through- 
put, and fringe modulation. There is one aspect 
of this figure (and also Figure [7]) that may ap- 
pear paradoxical, and so deserves some clarifica- 
tion: We would expect the fringe period (measured 
e.g. in microns) will be larger at the longer wave- 
lengths, which means that the fringes should be- 
come progressively more "stretched out" at larger 
X. However, it is clear that the fringes in Fig- 
ures |6] and [7] are actually more "bunched up" at 
large X. The reason for this can be found in 
the dispersion curve of Figure [H this curve shows 
that the same interval in X covers a much larger 
range in wavelength at the red end of the spec- 
trum than it does at the blue end. The results are 



^It could as well be perpendicular to the baseline, or at any 
other arbitrary orientation. 



clear from the labelling of the x axes in Figures [5] 
and [6l For instance, the 500 interval in X from 
< X < 0.5 mm covers a 50 nm wavelength in- 
terval from « 400 — 450 nm at the blue end (left), 
whereas the same interval from 1.5 < X < 2.0 
mm covers a wavelength interval that is w 6 times 
larger at the red end (right), from w 700 — 1000 
nm. It follows that in these figures the fringes 
will appear to bunch together at large X, even 
though the intrinsic fringe period is a factor of 
w 2.5 greater there. 

4.4. Delay scanning 

As described in Paper I, the fringe measure- 
ment process for SIM consists first of mechanically 
adjusting a coarse component of the delay in or- 
der to position the white fringe near to the tar- 
get. The internal delay S is then scanned over a 
small range, typically ±1A, and data recorded at 
the various wavelengths. This scanning of the de- 
lay leads to changes in the dispersed images owing 
to the varying modulation of the fringe function. 
Figure [7] shows the dispersed images of the "tar- 
get" and the "field star" in Figure [S] at 3 settings 
of the internal delay, from -500 to -1-500 nm. For 
this simulation, the target star is assumed to cor- 
respond to OPD = 0, so that the middle panel 
shows the SED of this star. At the extremes of 
the internal delay, slow changes in the spectrum 
of the target star are present, whereas the general 
appearance of the field star does not seem to vary 
much. However, the resolution in the focal plane 
(24 iim, see next section) is sufficiently fine that 
rapid variations can sometimes be seen from one 
wavelength channel to the next. 

4.5. Binning 

The final step in our model of the "real" focal 
plane is to account for the size of the CCD pixels, 
and for the integration in each CCD pixel during 
scanning of the internal delay. Both these steps 
amount to binning during the acquisition of the 
data. 

4-5.1. Pixellation in the focal plane 

The CCD in the focal plane will have 24 x 24 ^m 
pixels, which we model by binning the 1 /im pix- 
els of our idealized focal plane simulation. The 
camera optics are such that the scale in the focal 



7 



Wavelength (nm) 

534 




1.0 1.5 

Focol plane position (mm) 



2.0 



Fig. 6. — A portion of the image in the ideal focal plane of SIM including all effects discussed in this paper: 
diffraction, dispersion, throughput, fringe modulation, and confusion. The FOV includes two stars, a G-type 
"target" star at the center, and an M-type "field star" offset by 1.5" towards the upper left corner of the 
FOV. In this simulation, A = on the target star, so we see only the SED of this star (times the throughput, 
which falls off at either end of the spectrum). The (faint) widely-spaced dark bands in this spectrum are 
spectral features of this star. Above this G star is the spectrum of the field M star which is shifted « 1" 
vertically upwards and w 1" horizontally to the left of the target star and, since it is offset substantially 
from the white fringe, it is also modulated by fringes determined from the fan diagram. 
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Fig. 7. — Focal plane intensity distribution for three different values of the internal delay. Only the central 
profiles in X of the target and the field star have been plotted. As the scanning progresses, both the target 
and the field star show variations at differing rates, depending on their (baseline-projected) locations in the 
FOV and the value of the internal delay. The throughput and the dispersion have been included, but not 
the diffraction (which will reduce the fringe contrast, especially on the field star). 



plane will be 2" per CCD pixel of 24 ^m. Fig- 
ure [8] shows the mapping between the intensity in 



the region of the 3" field stop and the CCD de- 
tector plane. The idealized focal plane is binned 
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Fig. 8. — Correspondence between the ideal focal plane, the field stop, and the CCD pixels in the FOV. 
Left: The dashed circle shows the 3" FOV, covered with 2" CCD pixels. Right: The diffracted image is 
dispersed in the +X direction over 80 pixels. 



to 3 X 80 pixels, which is the effective sensitive 
area of the SIM CCD to be used for analysis of 
the fringes. The current default is to carry out 
a summation of the data in rows 1, 2, and 3 for 
each of the 80 channels, and to provide for a se- 
ries of on-board binning patterns of the data in 
the "channel" direction. A different approach of 
reading out each row separately may confer an ad- 
vantage in crowded fields; we will have more to say 
about this later in this paper. 

4-5.2. Integration of the fringe signal 

In addition to the pixcUation introduced in the 
previous section, the fringe signal of our ideal- 
ized focal plane is integrated over a finite delay 
range (and time) in order to simulate the opera- 
tion of the fringe detector. This integration oc- 
curs during scanning of the internal delay. The 
scanning involves moving a dither mirror back and 
forth with an amplitude of 1 /im on top of an ad- 
ditional fixed delay following a triangular-shaped 
wave with a nominal frequency of 250 H^, and si- 
multaneously recording eight or sixteen exposures 
of 250 IIS in each of all 80 channels. The total 
power Pq, the resultant fringe amplitude A, and 
the resultant fringe phase 4> are then determined 



*The integration time of the fringe detection system can be 
lengthened simply by decreasing the scanning frequency. 
This facilitates observations of faint targets with SIM. 



for each channel using a least squares method. In 
our simulation of this process we calculated ide- 
alized focal plane intensity distributions for ev- 
ery nanometer over the internal delay range of 
—0.5 < 5 < +0.5/im and binned them to obtain 
8 samples, each corresponding to a delay step size 
of exactly 125 nm. We then used t he on-board 
algorithm ( Milman fc Basingei 2002f) in order to 
calculate the fringe parameters. 

5. A comparison of the models 

In Paper I, we applied our simplified model to a 
number of SIM Key Projects and estimated the ex- 
pected single-measurement astrometric bias aris- 
ing from confusion. In that paper, we also tab- 
ulated the differences with the results using the 
improved model of the present paper for several 
representative cases. Table [T] reproduces those re- 
sults here for convenience. We have already noted 
that the differences are not significant, so that a 
good estimate of the presence and magnitude of 
the confusion bias can be obtained without the 
additional complexity of the improved instrument 
model. In this section we discuss the comparison 
between the two models in more detail. 

In Paper I we found that the magnitude of the 
confusion bias depends on several factors: the rel- 
ative brightnesses of the field stars with respect 
to the target star (including any vignetting by the 
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entrance aperture); the SEDs of the target and 
field stars; and the separation of the field stars 
from the target star as projected along the inter- 
ferometer baseline. Since these parameters have 
not changed with the advent of the improved in- 
strument model, and since the calculations have 
been carried out without the addition of any noise 
sources, the results ought to be identical. How- 
ever, closer inspection of Table [T] shows small but 
systematic differences, as follows: 

• The total confusion biases we give in Table 
[1] for the sample fields computed using the 
improved model (column 4) are all smaller 
than those of the simple model (column 3), 
except for case 5 (where the difference is 
zero) , and case 7 (where the improved model 
is worse by 2 ^as). 

• The largest (absolute) differences (column 5) 
occur for cases 2, 4, and 1, in that order. 
Case 2 involves a blue Z—0 quasar and a 
very red field star of type M. Case 4 has all 
the same instrument parameters as case 2, 
and also involves a very red M field star, but 
the quasar is now at Z= 2, i.e. it is less blue. 
Case 1 is identical to case 2, but now the 
field star is a less-red type A. 

To summarize, the improved model generally 
predicts smaller astrometric biases than does the 
simple model, and the largest differences between 
the models appear to occur for cases involving the 
reddest field stars. We now look for the explana- 
tion of these effects in the differences between the 
models. 

5.1. Differences in the models 

The improved model uses the same phasor de- 
scription as the simple model, and consideration of 
this description provides a basis for understanding 
the origin of the differences. Figure 6 of Paper I 
shows an example of a confused field. In that fig- 
ure, the target phasor is shown as the thick black 
arrow pointing along the -fX (real) axis, and 10 
additional faint confusing stars are scattered over 
the field. The resulting phasor no longer points 
along the -l-X axis, but has small errors both in 
amplitude and phase compared to the target. The 
phase error causes the astrometric confusion bias. 
Now suppose we have just one confusing field star; 



the only way to reduce the amount of the bias (j) 
in that case is to reduce the amplitude of the pha- 
sor representing that field star. So we must look 
for any differences in the improved model which 
might result in reducing the fringe amplitude of 
the field star. 

The first, and perhaps most obvious, way of 
reducing the contribution of the field star is to re- 
duce the size of the pixels. The improved model ef- 
fectively has 2" X 6" square pixels, which is a factor 
1.7 larger than the 3" round pixels of the simple 
model, so this unfortunately goes the wrong way. 
However, with the exception of case 6, the angular 
separation between the target and the field star is 
a small fraction (« 50/300 = 0.17) of the collector 
PSF, so this effect is unlikely to contribute much. 
The exception is case 6 in Table [T] (angular sepa- 
ration 1.5"), but this has the faintest field star (at 
3 mag below the target) of all the cases in Table, 
so its bias effect was negligible to start witt|f|. 

The next difference between the models to con- 
sider is that part of the "binning" which involves 
integration of the fringe signal during the de- 
lay scanning. However, this refinement is un- 
likely to have any effect at all, for the follow- 
ing reason: The simple algorithm modeled the 
fringe as a densely-sampled sinusoid, while the 
improved model uses a much coarser binning of 
the fringe into eight segments. However, the al- 
gorithm used to compute the fringe parameters is 
the same for both models; it is a standard ap- 
proach which explicitly accounts for the time bin- 
ning using the fact t hat the signal is a s i nusoi d 
of known period (e.g. Milman fc Basingen ([2002), 
Catanzarite fc MilmanI (|2002l )). This refinement 



of the improved model should therefore have no 
effect. 

Finally, there is the binning of the dispersed 
fringes which occurs owing to the pixellation of 
the CCD in the focal plane. This remaining effect 
appears capable of accounting for the differences 
we have observed. The situation can perhaps be 
best understood by referring to Figure [T] This fig- 
ure shows that the period of the dispersed fringes 
is very long for the target star, but much shorter 
for field stars. The example shown in this Figure 



^We will return later to the suggestion that reducing the 
(effective)pixel size on the CCD detector can be a way of 
reducing the confusion bias for some special cases. 
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is an extreme case of course, with a projected sep- 
aration of « 1". After pixellation by the CCD, the 
amphtude of the target star fringe which wiU be 
recorded in any given pixel (channel) as the delay 
is scanned will hardly be affected, because the pe- 
riod of the dispersed fringes is many dozens of pix- 
els (of order 80 if the target is at the white fringe). 
But the periods of field stars are shorter, and the 
extreme example shown in this figure has a period 
of about 1 pixel near the middle of the CCD. It's 
clear that the fringe amplitude measured on this 
field star as the delay is scanned will be greatly 
reduced, indeed, essentially zero in this example. 
Field stars not as distant from the target will be 
less affected, but the amplitudes of their phasors 
will nevertheless be smaller, leading to a reduc- 
tion in the confusion bias they contribute. Note 
that this effect will never lead to an increase in the 
amount of the confusion bias, only a decrease, as 
is observed in Table [TJ To conclude the argument, 
we note that the period of the dispersed fringes 
on the CCD (cf. Figure [6]) is larger at the blue 
end of the spectrum (about 2 pixels in this exam- 
ple) than it is at the red end (about 0.8 pixels). 
The smoothing effect of pixellation on the CCD 
will therefore be most severe for the reddest field 
stars, as is also observed in Table [1] We conclude 
that it is the pixellation by the CCD in the focal 
plane which is responsible for the reduction in the 
level of confusion bias recorded with the improved 
instrument model described in this paper. 

6. Discussion 

A number of issues have been uncovered dur- 
ing the course of this work which deserve further 
scrutiny, but for which there is insufficient space 
in this paper to deal with properly. We briefly 
mention several of these remaining issues. 

6.1. Separate row readout 

As explained earlier, (cf. Figure [8]), there will 
be three rows of 80 pixels on the detector devoted 
to collecting the dispersed fringe photons. The 
nominal SIM design is to sum the charge from all 
3 rows on board to form a single row of 80 pix- 
els; the fringe parameters will be estimated from 
these data as the internal delay is scanned. Here 
we show that, in crowded fields, it may be useful to 
read the rows out separately, and to estimate three 



sets of fringe parameters, one from each row. For 
example, consider a field consisting of a binary in 
which one of the stars (which we designate as the 
"target") is in the middle pixel, and the other star 
(the "field" star) is in the upper pixel. The field 
star in this simulation is fainter by 2 magnitudes, 
and located at a radial distance of 1.4" at a posi- 
tion angle of 2°. The projected separation (along 
the interferometer baseline) is about 50 mas. Fig- 
ure ini shows the signal spectrum on the detector 
with and without the field star. This suggests that 
the astrometric bias arising from confusion would 
be significantly reduced if only the central row of 
pixels is used to compute the fringe parameters. 
A comparison of the phase spectrum computed 
for each case confirms this suggestion, as shown 
in Figure [10] (See Paper I for the definitions and 
a discussion of the phase and delay spectra) . The 
phase spectrum computed using only the central 
row of pixels shows the least perturbation, and the 
resulting bias in the astrometric delay caused by 
this field star is small (0.17 nm). However, using 
all 3 rows results in a significant delay bias which 
is almost 4 times larger. Clearly, in this case the 
best strategy would be to use only the central row 
of pixels. In fact, looking for consistency in the 
fringe parameters calculated from each row sepa- 
rately may be a useful indicator of the presence 
of confusion, especially if the details of the distri- 
bution of field stars is not known ahead of time. 
This subject clearly warrants a more detailed dis- 
cussion; it is possible that the advantages of read- 
ing out the data for each row separately may out- 
weigh any signal-to-noise penalty. In Paper I we 
concluded that modeling and removing the astro- 
metric bias arising from confusion in SIM mea- 
surements is not likely to be feasible owing to the 
lack of sufficiently-detailed SEDs and accurate po- 
sitions of all the relevant stars in the FOV. How- 
ever, whatever level of bias may be present in the 
data, it is quite possible that this bias can be fur- 
ther reduced by using the additional flexibility of- 
fered by separately analyzing the 3 rows of the 
CCD detector. 

6.2. Self-confusion 

A description of how SIM can in principle 'con- 
fuse itself was given in §4.21 The basic prob- 
lem here is that, even when there is only one star 
within the FOV, the effective wavelengths of the 
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Table 1 

Comparison of the single-measurement astrometric bias o btained with the improved 

INSTRUMENT MODEL TO THE RESULTS OF THE SIMPLIFIED MODEL OF ISRIDHARAN fc AlLEnI (j2007h . 



Case 


Model 


Simple 


Improved 


1 Difference! 


Number 


parameters 


(/ias) 


{fia.s) 


(flSLS) 


1 


Q(z=0), AlV, 2, 50, 5, 90 


442 


431 


11 


2 


Q(z=0), M6V, 2, 50, 5, 90 


1072 


1036 


36 


3 


Q(z=2), AlV, 2, 50, 5, 90 


224 


215 


9 


4 


Q(z=2), M6V, 2, 50, 5, 90 


614 


592 


22 


5 


Q(z=2), M6V, 2, 50, 5, 10 


-46 


-46 





6 


AlV, BIV, 3, 1500, 10, 90 


-1 





1 


7 


AlV, M6V, 2, 50, 90, 90 


-59 


-61 


2 


8 


AlV, M6V, 2, 25, 90, 90 


-53 


-46 


7 



Note. — In column 2, Q(z=2), AlV, 2, 50, 5, 90 refers to a model with a 
redshift 2 target quasar and an AlV field star with Am = 2, located 50 mas 
distant in PA = 5°. The baseline orientation is 90°. Cases 1-5 are taken from 
the Quasar Frame Tie key project (cf. Paper I). Cases 6-8 are taken from a 
selection of binary models. 
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Fig. 9. — T(MR) (solid line); signal spectrum when only a single target star is present at the delay center 
in the FOV, and only the middle row of pixels is used. T-fFS(MR); as previously, but with the field star 
present. T-fFS(AR); as previously, but adding the signal in all 3 rows. Note how the contribution from the 
field star has grown in the last plot, dominating the signal around channels 30 and 55. 
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Fig. 10. — Phase spectrum obtained separately from each of the three rows of the CCD detector for the 
binary target described in the text, and finaUy from the summed signaL In aU 4 cases the bias in the 
astrometric delay caused by the presence of the field star is included. The bias obtained from adding the 
intensities of all the three rows (0.64 nm) is higher than that determined using only the middle row(0.17 nm). 



channels will be affected by photons diffracted into 
that channel from neighboring channels. This ef- 
fect will clearly depend on the spectrum of the 
target star, including e.g. absorption features. We 
find that there could be a maximum of 0.5% 
change in the effective wavelength of a given chan- 
nel both for A and M type stars. Accordingly, we 
estimated the fringe phase first with the actual 
wavelengths of each channel, and then with a ran- 
dom perturbation of those wavelengths by 0.5%. 
The differences in the astrometric delay are 0.2 nm 
for an A type star and 0.35 nm for an M type star. 
The corresponding single-measurement bias in the 
position measurement is « 8^as, which is likely to 
be negligible in most cases. 

The question arises as to just how the effective 
wavelength of each channel will be calibrated in 
orbit. If a bright star is used for this purpose, 
the effects of 'self-confusion' described above will 
be included in the calibration. However, it is also 
clear that the results of such a calibration will be 
dependent on the SED of the calibrating star, so 
this choice will have to be made carefully. 



There is a variant on this problem when the 
FOV includes field stars as well as the target star. 
In this case it is clear that the effective wavelength 
of the sum of all the photons arriving at any po- 
sition in the focal plane will depend on the SEDs 
of all the field stars as well as that of the tar- 
get star. Our simulations indicate that this prob- 
lem is slightly more severe than the 'self-confusion' 
described above, and in general the change in 
wavelength arising from the presence of multiple 
sources is greater in the longer wavelength chan- 
nels. We have estimated this effect in a few test 
cases and find the biases to be in the range of 
1-10 nm (20-200 /las) at short wavelengths, and 
slightly more at long wavelengths. The pixellation 
of the focal plane tends to smooth out this change 
in effective wavelength to some extent. Again, one 
way at least to identify the presence of a possible 
problem is to estimate the fringe parameters in 
each row separately; if they are inconsistent, it 
may help to ignore the red end of the phase spec- 
trum. 



Milman et al.l (|2002l ) have modeled the path- 
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delay error for SIM arising from wavelength errors 
and shown it to be small. However, their analy- 
sis did not include the diffraction effects we have 
described above. 

6.3. Loose ends 

Finally, we note that even our "improved" in- 
strument model contains numerous idealizations 
which deserve further scrutiny. For example, our 
model for the pixellation of the focal plane as- 
sumes there is no space between the CCD pixels, 
and the response across these pixels is uniform. 
Similarly, variations of parameters across the FOV 
such as the throughput and the dispersion have 
been ignored. 

7. Conclusions 

We have developed an improved instrument 
model for SIM which includes a number of impor- 
tant effects that were not considered in our first 
paper. The additional features included here pro- 
vide a much more detailed understanding of just 
how SIM works. We note that this significantly- 
more-complicated model does not lead to any ma- 
jor changes in our initial estimates of the astromet- 
ric bias arising from extraneous field stars in the 
field of view (cf. Paper I), and the small differences 
which have appeared for the sample fields we have 
calculated can be explained as a consequence of 
smoothing by the pixellation of the images on the 
CCD in the focal plane of the camera. Besides the 
obvious pedagogical value of the improved model, 
it is likely to be a useful point of departure for 
future investigations of even more subtle instru- 
mental effects which may be present in real SIM 
data. 

We are grateful to our colleagues in the SIM Sci- 
ence and Engineering Teams for discussions about 
the potential effects of confusion on SIM astrom- 
etry. Special thanks go to Mike Shao for sharing 
his detailed knowledge about the design of SIM, 
and to him and Steve Unwin for their encourage- 
ment to include the effects of the measurement 
process discussed here into our modeling of confu- 
sion as described initially in Paper I. We have es- 
pecially appreciated the constructive criticisms of 
the referee on an earlier version of this paper. This 
work was funded by the Jet Propulsion Labora- 



tory under contract ^1268384, and carried out at 
the Space Telescope Science Institute. We thank 
the Institute director for providing additional fi- 
nancial support to permit the timely completion 
of this project and publication of the results. 

For queries on the IDL code developed for these 
studies please contact the authors. 
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A. Modeling the SIM focal plane 

Wc have chosen to model the ideal focal plane with a X,Y rectangular grid of 2112 x 192 elements, based 
on the following considerations: 

1. The spatial sampling at the focal plane was chosen for convenience to be 1 /im. The CCD camera will 
have 24 x 24 ij,m. square pixels, each of angular size 2" x 2" on the sky. The linear size of each array 
element in our model of the focal plane therefore corresponds to an angular size of 1/12". 

2. In order to include the diffracted light from a star located as much as 3" from the center of the FOV, its 
point spread function was simulated with a separate array of grid size 120 x 120 elements, corresponding 
to a total angular size of 10" x 10". 

3. We have required that the light from all the stars within a circle of diameter 6" centered on the FOV 
be included in the simulation. This implies that the minimum size of the simulated FOV should be 
6" X 6". Accounting for the size of the point spread function, the required simulation area grows to 
16" X 16". The simulated focal plane must therefore be 16 x 12 = 192 elements in size along the 

direction normal to that of the dispersion (the Y-direction) . 

4. The design of the camera was such as to provide 80 wavelength channels, i.e. 80 CCD pixels along the 
direction of the dispersion. This corresponds to a linear size of 1920 microns. 

5. The spectra of stars not at the center of the FOV will be shifted along the direction of dispersion. This 
may require as much as an additional ±36 elements in order to include all the relevant stars within the 
circle of diameter 6". Thus, the required number of array elements along the direction of dispersion in 
+X is 1920 + 2 X 36 + 120 = 2112. 
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